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ABSTRACT 


Assuming  that  Saturn  has  trapped  radiation  belts  similar  to 
those  of  Earth  and  Jupiter,  it  is  shown  that  observations  of  the 
radial  dependence  of  the  intensities,  energy  spectra,  electron- to- 
proton  intensity  ratio,  and  pitch  angle  distributions  of  energetic 
charged  particles  trapped  outside  of  Ring  A can  potentially  pro- 
vide information  -fa^T  on  the  existence  of  Ring  E,  {V f on  the  size 
distribut  on  of  the  particulate  matter  therein,  and  (c)'  on  the 
magnitude  of  the  radial  diffusion  coefficient  for  energetic 
particles  A parametric  study  of  these  possibilities  is  specialized 
to  the  characteristics  of  the  University  of  Iowa  detectors  on 
Pioneer  11,  which  will  make  a close  encounter  with  Saturn  in 
August-September  1979- 
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INTRODUCTION 


During  the  1966  edge-on  presentation  of  the  ring  plane 
of  Saturn  to  Earth,  Feibelman  [1967]  obtained  photographic 
evidence  for  a very  faint  ring  external  to  the  well  known  visible 
rings.  The  evidence  for  this  external  ring  (first  called  Ring  D' 
and  now,  by  I.A.U.  convention,  Ring  E)  las  been  reexamined  by 
Smith  et  al.  [1975]  and  Smith  [1978].  Eased  on  the  latter 
reference,  it  appears  that  the  normal  optical  opacity  T]  of 
Ring  E lies  in  the  range  0.5  X 1CT7  to  3 X 10"6,  if  full 
credence  is  given  to  Feibelman' s photog -aph  of  15  November  1966, 
though  there  is  no  observational  knowledge  whatever  on  the  size 
distribution  of  the  particulate  matter  therein.  (Note  that  T] 
is  of  the  order  of  unity  for  Rings  A and  B. ) 

The  question  of  the  existence  of  a Saturnian  magnetic 
field  and  of  the  nature  of  a Saturnian  magnetosphere  has  been  dis- 
cussed by  Van  Allen  [1977]  and  Siscoe  [1978a,  b].  For  the  present 
analysis  we  assume  that  Saturn  has  a magnetosphere  which,  at 
least  near  the  planet,  is  similar  to  that  of  Jupiter;  that  is, 
we  assume  a dipolar  magnetic  field  configuration  and  a significant 
population  of  trapped  particles  for  which  the  source  lies  outside 
the  inner  magnetosphere  and  for  which  the  dominant  transport 
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mechanism  is  radial  diffusion  violating  the  third  adiabatic 
invariant  while  the  first  and  second  are  conserved.  Such  a picture 
has  enjoyed  considerable  success  in  providing  an  understanding  of 
the  inner  magnetosphere  of  Jupiter  [Schulz,  1979]* 

The  ring  material  within  Saturn's  magnetosphere  is  a sink 
for  trapped  radiation,  with  the  lifetime  against  absorption  by  the 
ring  being  a function  of  charged  particle  species  and  energy  as 
well  as  a function  of  the  ring  properties.  As  discussed  by 
Van  Allen  [1977]>  Kings  A and  B are  sufficiently  opaque  that  the 
lifetime  against  absorption  by  them  is  of  the  order  of  a latitudinal 
bounce  period,  so  that  they  may  be  treated  as  perfect  barriers  to 
the  inward  diffusion  of  charged  particles.  Ring  E,  however,  may 
be  sufficiently  transparent  that  a significant  population  of 
charged  particles  is  able  to  survive  long  enough  to  diffuse  into 
the  outer  edge  of  Ring  A,  at  which  complete  absorption  occurs. 

In  this  paper  we  study  parametrically  the  absorption  of  charged 
particles  by  Ring  E in  order  to  investigate  its  effect  on  charged 
particle  intensities,  spectra,  and  pitch  angle  distributions, 
all  of  which  are  potentially  observable  oy  the  energetic  charged 
particle  instruments  on  Pioneer  11. 
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PROCEDURE 

Our  procedure  is  to  assume  that  particles  are  injected  with 

a specified  energy  spectrum  and  pitch  angle  distribution  at  some 

outer  boundary,  which  for  convenience  we  take  to  be  at  the  dipole 

L- value  of  5-  (All  radial  distances  are  given  in  units  of 

Saturn's  equatorial  radius,  60,000  km.)  The  particles  are  then 

assumed  to  diffuse  radially  inward  until  they  are  absorbed  at  the 

outer  edge  of  Ring  A at  L = 2.27.  The  diffusion  coefficient 

governing  the  radial  motion  is  taken  to  be  of  the  form  D = D Ln. 

o 

In  analogy  to  Jupiter's  magnetosphere,  we  adopt  the  value  n = 3, 

which  is  appropriate  for  diffusion  driven  by  ionospheric  winds, 

the  driving  mechanism  which  appears  to  dominate  in  Jupiter's  inner 

magnetosphere  [Schulz,  1979]-  Absorption  of  charged  particles 

by  Ring  E material  is  assumed  to  be  the  only  loss  process,  with 

lifetime  t Lm. 
o 

The  differential  unidirectional  Intensity  at  energy  E, 
equatorial  pitch  angle  a,  and  dipole  L- value  L is 

j(E,a,L)  = pa  f(n,Ja,L)  (1) 

where  p is  the  momentum  of  a particle  with  energy  E,  n(E,a,L)  is 
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the  first  adiabatic  invariant,  J2(E,a,L^  is  the  second  adiabatic 
invariant,  and  the  phase  space  density  f is  the  solution  to  the 
radial  diffusion  equation 


La 


A. 

bh 


"l 


J 


(2) 


with  boundary  conditions 


f(L  = 2. 27)  = 0.0 


(3a) 


f(l  - Lq  - 5) 


iol*o’ao’Lo)W  ‘ 


(3b) 


E , p , and  a are  related  to  E,  p,  and  a by  conservation  of  the 
o o o 

first  and  second  invariants: 

E0  “ Eo(n,Ja,L) 

Qq  * <*0(n> Ja>L)  • 

The  boundary  condition  (5b)  is  equivalent  to  specifying  the  differen 
tial  unidirectional  intensity  at  L * 5 as  a function  of  ctQ  and  Eq. 
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LIFETIMES  AGAINST  RING  PARTICLE  ABSORPTION 

For  ring  particle  absorption,  the  quantities  tq  and  m 
in  equation  (2)  can  be  determined  as  functions  of  a,  the  ring 
particle  radius  r , and  the  optical  opacity  of  the  ring  71  in  the 
following  way. 

The  average  thickness  of  a spherical  ring  particle  with  radius 

r is  4r  / 3 . The  quantities  t and  m depend  on  whether  the  range  R 
o o o 

of  the  charged  particle  in  the  ring  material  (which  we  assume  is 
water  ice)  is  greater  than  or  less  than  this  average  thickness. 

R > 4ro/3 

The  average  time  for  a charged  particle  to  make  a single 
encounter  with  a ring  particle  is 

Ti  = TQ  cos  a/27]  (4) 

where  Tn  is  the  latitudinal  bounce  period  of  the  charged  particle. 
Since  the  average  thickness  of  ring  material  traversed  in  a single 
encounter  is  4r J3,  the  average  thickness  traversed  per  unit  time 
is  871  r ^TgCosa  . If  we  define  the  lifetime  against  absorption  as 
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the  time  it  takes  a particle  traversing  ring  material  at  this  rate 
to  traverse  a distance  equal  to  its  range,  we  have 


T a 3Tfi  R cos  a / 8T1  rQ  . 


(5) 


To  find  tq  and  m,  we  need  to  express  Tg  and  R as  functions  of  p, 
and  L.  For  this  purpose  we  write  the  particle  energy  as 


E 


(j3Q/L3sin2  a 
( an^c^B^/L3  s in2a ) ^ 


where  Bq  is  the  equatorial  surface  fieli  strength  and  mQc2  is  the 
rest  energy  of  the  particle,  (in  this  and  subsequent  relations, 
the  upper  expression  is  for  non-relativistic  particles,  and  the 
lower  one  is  for  ultra- relativistic  particles.) 


TB 


4r  Lf(sina)  x ( 


(m0L3  sin2  a /2p  Bq  )* 

L/c 


(6) 
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where  f(sina)=  1.3802  -0.3199  [sina  + /sina  ’]  [Schulz  and 
Lanzerotti,  19?U].  In  addition,  we  can  approximate  the  particle 
range  as  a power  law  in  energy  (which  is  valid  over  fairly  wide 
ranges  of  energy),  i.e. 


R « A Ea 

r 

(p,B  Jl?  sina  a)a 

« A ( 

j (2moc2liBo/L3  sinaa)a/a 


(7) 


Combining  (6)  and  (7)  we  find 


m 


and 

T 

< 

where 

G 

H 

g 

h 


f G(rQT|)  «(a)  (^B0)a"^ 
Y H(roTl)  h(a)  (^Bo)a/2 


3ARg  (imoca)^/2c  tqT\ 
3ARsJ(2moca)a/s/2c  rQ  T) 
cosa  sin1_aa  a f(sinCt) 
cosa  sin  aa  f(sina) 


(8) 

(9) 
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PARAMETERS  WHICH  GOVERN  THE  SOLUTION 

The  solution  to  equation  (2)  is  given  by  equation  (8)  in 

Thomsen  el.  al.  [1977].  As  discussed  in  that  paper,  the  quantity 

which  governs  the  nature  of  the  solution  is  D^t^.  This  is 

because  D - 1/ (diffusion  time)  so  that  DqTo  is  roughly  proportional 

to  the  ratio  of  the  loss  time  to  the  diffusion  time.  If  DqTo  is 

very  large,  the  particles  suffer  few  losses;  if  ^0T0  small, 

the  losses  are  severe.  From  the  above  expressions  for  tq,  we  see 

that  the  quantity  which  governs  the  solutions  for  the  ring  sweep-up 

problem  is  Do/rQT)  if  R > 4ro/3  and  is  just  Dq/T1  if  R < 4rQ/3. 

In  Figure  1 we  explore  qualitatively  the  dependence  of  the 

solution  on  the  ring  particle  radius,  the  ring  opacity,  and  Dq  for 

the  various  species  and  energies  of  charged  particles  which  are 

observable  by  the  University  of  Iowa  detectors  on  Pioneer  11. 

(For  a description  of  the  instrument,  see  Van  Allen  et  al. , 1975- ) 

The  curves  plotted  in  Figure  1 show  the  value  of  which  results 

in  Dt  = D t Lm+n  = 1 at  L = 3.5  for  each  value  of  the  ring  particle 
o o 

radius.  Since  the  condition  Dt  = 1 means  that  the  diffusion  time 
is  of  the  same  order  as  the  lifetime  against  loss,  the  lines  shown 

i : 

in  Figure  1 roughly  separate  the  range  of  parameters  for  which 
losses  are  severe  from  the  range  for  which  losses  are  negligible. 

I 


4 
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Thus,  for  example,  if  the  ring  particle  radius  is  1 cm  and  D /T) 

has  the  value  10"1,  then  we  would  expect  protons  with  an  energy 

of  0.6l  MeV  to  suffer  very  few  losses,  while  electrons  with  an 

energy  of  O.56  MeV  should  be  significantly  depleted.  In  general, 

particles  whose  curve  in  Figure  1 lies  highest  for  a given  value 

of  rQ  will  be  the  most  strongly  depleted.  Thus  below  rQ  « 0.4  cm 

protons  with  an  energy  of  0.6l  MeV  are  more  strongly  depleted 

than  electrons  with  an  energy  of  21  MeV,  whereas  above  that  energy 

the  reverse  is  true.  Therefore,  if  Pioneer  11  observes  that 

21  MeV  electrons  are  more  strongly  depleted  than  are  0.6l  MeV 

protons,  we  know  that  is  probably  greater  than  0.4  cm. 

Examination  of  Figure  1 shows  that  the  crossing  points  for  the 

various  curves  occur  at  different  values  of  r so  that  by  inter- 

o 

comparing  observations  made  by  the  several  detectors,  one  may  be 
able  to  bracket  the  value  of  r . 


SAMPLE  RESULTS 


In  Figures  2-4  we  investigate  the  effect  of  ring  particle 
absorption  on  the  integral  intensities  observable  by  Pioneer  11 
detectors.  To  obtain  these  figures  we  have  assumed  that  the  uni- 
directional differential  intensity  at  Lq  = 5 is  isotropic  and  is 
a power  law  in  energy:  jQ(Eo,  aQ>  4q)  = K E“2.  We  have  then 
followed  the  procedure  described  above  to  determine  j(E,  a,  L) 
over  a range  of  energies  at  each  value  of  L.  The  integral 
intensity  corresponding  to  each  detector  is  found  by  integrating 
j(E,  a,  L)  over  the  energy  range  of  the  detector 

Figure  2 shows  the  radial  profiles  of  unidirectional  integral 
intensities  that  would  be  observable  at  a local  pitch  angle  of  85° 
on  the  inbound  pass  of  the  Pioneer  11  Saturn  encounter  if  the 
Saturnian  magnetic  field  were  a dipole  aligned  with  the  rotation 
axis  of  the  planet.  The  three  profiles  shown  for  each  detector 
correspond  to  the  loss-free  case  and  to  the  two  sets  of  the  ring 
particle  parameters  marked  with  an  x in  Figure  1.  Note  how  the 
extent  of  the  deviation  from  the  loss- free  profile  depends  on  how 
far  below  the  curves  of  Figure  1 the  adopted  parameter  set  lies. 

Figure  3 illustrates  the  effect  which  ring  absorption  has  on 
the  energy  spectrum  and  proton/electron  intensity  ratio  for  the  three 
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cases  illustrated  in  Figure  2.  Shown  in  Figure  3 are  the  ratios  of 
intensities  that  would  be  seen  by  three  of  our  detectors  to  the 
intensity  that  would  be  seen  by  the  detector  with  an  electron  energy 
threshold  of  40  keV.  As  expected  from  figure  1,  the  effect  of  the 
ring  absorption  for  the  two  cases  shown  is  to  deplete  low-energy 
electrons  more  than  high-energy  electrors  or  low-energy  protons. 

Figure  4 illustrates  the  effect  which  ring  absorption  has  on 
the  pitch  angle  distribution  of  the  particles,  again  for  the  three 
cases  illustrated  in  Figure  2.  In  this  figure  are  shown  the  ratios 
of  the  intensities  at  a local  pitch  angle  of  85°  to  those  at  a local 
pitch  angle  of  45°.  It  is  evident  that  there  is  a strong  pitch  angle 
dependence  to  the  ring  absorption  losses,  with  large  pitch  angle 
particles  suffering  more  severe  losses.  The  result  is  a strong 
tendency  for  the  pitch  angle  distributions  to  become  more  dumbbell 
in  nature  (i.e.,  peaked  near  0°  and  l8or')  at  smaller  radial  dis- 
tances. This  is  opposite  to  the  effect  of  loss-free  radial 
diffusion,  which  is  for  the  distributions  to  become  more  pancake 
in  nature  (i.e.,  peaked  near  90°).  This  effect  on  the  pitch 
angle  distributions  is  probably  the  best  diagnostic  for  the 
presence  of  ring  particle  sweep-up  as  a significant  loss-process. 
(Note,  however,  that  the  effect  of  a tilt  between  the  dipole  and 
rotational  axes  is  to  increase  the  lifetime  of  particles  mirroring 
below  the  maximum  magnetic  latitude  excursion  of  the  ring  [see,  e.g., 
Mead  and  Hess,  1973]  and  thereby  reduce  the  strength  of  the  effect 
on  the  pitch  angle  distributions. ) 
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SUMMARY  AND  CONCLUSION 

To  summarize  the  results  of  the  above  analysis: 

(1)  The  presence  of  a tendency  Tor  pitch  angle  distribu- 
tions (especially  of  low  energy  electrons)  to  become  more  dumbbell 
at  smaller  radial  distances  is  probably  the  clearest  potential 
evidence  that  Ring  E absorption  losses  are  significant. 

(2)  Other  effects  of  diagnostic  value  are  the  inward 

radial  dependence  of  the  electron  energy  spectrum  and  the  intensity 
ratio  of  protons  0.6l  < < 3.41  MeV  to  electrons  Eg>  40  keV. 

As  shown  :.n  the  several  figures,  there  is  a significant  possibility 
for  determining  the  approximate  value  of  the  mean  radius  rQ  of 
ring  particles. 

(3)  If  rQ  can  be  determined,  the  degree  of  the  deviation 

of  the  observed  profiles  from  the  loss-free  case  can  be  used  to 

infer  an  approximate  value  of  Dq/T].  Since  7)  is  found  observa- 

tionally  to  lie  in  the  range  0.5  X 10-7  < Tl<  2.3  x 10-6 

[Smith,  1978],  one  may  then  be  able  to  place  limits  on  the  value 

of  D . 
o 

In  conclusion,  we  have  demonstrated  that  observations  of 
the  radial  dependence  of  the  intensities,  energy  spectra,  electron- 
to-proton  intensity  ratio,  and  pitch  angle  distributions  of 
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energetic  charged  particles  trapped  in  a Saturnian  magnetic  field 
outside  of  Ring  A can  potentially  provide  information  on  the 
existence  of  Ring  E,  on  the  size  distribution  of  the  particulate 
matter  therein  and  on  the  magnitude  of  the  radial  diffusion 
coefficient  for  energetic  particles. 
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FIGURE  CAPTIONS 

Figure  1.  Values  of  Do/T)  for  which  the  diffusion  time  (~  l/D)  is 

roughly  equal  to  the  lifetime  against  ring-particle  loss  (t) 
at  L = 3.5  for  the  various  species  and  energies  observable 
by  the  Pioneer  11  University  of  Iowa  detector  complement. 

For  Ee  = 0. 56  MeV  neither  the  non-relativlstic  limit  nor  the 
ultra-relativistic  limit  for  the  lifetime  (eqs.  8-11)  is 
appropriate,  so  the  lifetime  (and  the  value  of  D plotted 
here)  has  been  calculated  using  both  limits,  and  the  area 
between  the  two  results  has  been  shaded.  If  lies  much 

above  the  curve  for  particles  of  a given  species  and  energy, 
those  particles  will  suffer  very  few  losses  from  ring  absorp- 
tion. If  D^TJ  lies  much  below  the  curve,  they  will  suffer 
significant  losses.  Points  marked  by  x correspond  to  param- 
eter sets  investigated  in  subsequent  figures,  and  the  shaded 
region  shows  the  predicted  ranged  of  D^T]  obtained  by  combining 
the  estimated  opacity  of  Ring  E [Smith,  1978]  and  the  radial 
diffusion  coefficient  scaled  from  its  range  of  values  at 
JUplter  [e.g. , Thomsen  et  al. , 1977]  according  to  the  scaling 
relation  described  by  Siscoe  [1978a]. 


Figure  2.  Radial  profiles  of  unidirectional  (local  pitch  angle 

* 8f°)  integral  intensities  observable  by  University  of  Iowa 
detectors  as  predicted  for  loss -free  radial  diffusion  and  for 
radial  diffusion  with  Ring  E absorption  for  the  two  sets  of 
ring  parameters  marked  by  an  x on  Figure  1. 

Figure  3.  Ratios  of  unidirectional  integral  intensities  observable 
by  various  University  of  Iowa  detectors  as  predicted  for  ring 
absorption  assuming  two  different  sets  of  Ring  E parameters. 

Figure  4.  Ratios  of  predicted  observable  unidirectional  integral 

intensities  at  pitch  angles  of  85 5 and  45*  for  two  different 


sets  of  Ring  E parameters. 


B-G79-354 

D0A)  for  which  Dr  = I at  L=  3.5  ( 0^=85°) 


NO  SIGNIFICANT 
LOSSES 


E0  = 0.04  MeV 


E_-  21  MeV 


SIGNIFICANT 

LOSSES 


igure 


